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Hybrid photonic nanostructures allow to engineer novel interesting states of light. One recent 
example of that are topological photonic crystals where nontrivial Berry phase of photonic band 
structure gives rise to topologically protected unidirectionally-propagating (chiral) edge states of 
photons. Here we demonstrate that by coupling an array of emitters to the chiral photonic edge state 
one can create collective strongly - correlated states of photons in a highly controllable way. These 
correlated states are topologically protected and have a number of remarkable universal properties. 
Namely, an outcome of scattering does not depend on the position of emitters and is given only 
by the universal numbers, zeroes of Laguerre polynomials; correlation functions demonstrate clear 
even-odd effect and the result of scattering off a system is robust with respect to inhomogeneous 
broadening. 



INTRODUCTION 

Light-matter interaction is around us in nature and 
has played a tremendous role in the development of cur- 
rent technology. Until recent decades it was sufficient to 
deal with this interaction on average, with many photons 
and many atoms involved. However, increasing minia- 
turization of basic constituents towards nano scale is a 
common trend of modern technology. Downscaling to 
single- atom and/or single-photon levels promotes some 
traditionally classical research areas into the quantum 
realm PQ, [2] [3], [I]- Therefore control over quantum pro- 
cesses of light-matter interaction will eventually be a vi- 
tal ingredient of emerging quantum devices and is as well 
important for the rapid developments in several related 
fields, including communication, signal processing, ultra- 
fast optics, optomechanical cooling, imaging and spec- 
troscopy and, of course, quantum information. Efficient 
manipulation and control however requires a relatively 
strong interactions on a level of single atom, single pho- 
ton and/or a single electron levels [5], [5J, [7], [E]> [5], [ID]. 
[IT] . This represent a significant challenge since a typical 
interaction scale for individual particles in our everyday 
life is given by the QED coupling constant a ps 1/137. 
Two possible ways to overcome this natural limitation 
and to increase the effects of correlations is either to use 
artificial materials and devices or to engage many-body 
interacting effects which would eventually lead to desired 
non-linearities in systems with reduced dimensionality. 

Recent tremendous experimental progress in fabricat- 
ing few-photon sources coupled to ID transmission lines, 
[12], [H], [H], PS, [El, [H]> PS opens a prospect for 
creating and manipulating strongly-correlated states of 
photons. This also triggered significant number of the- 
oretical studies [T5],[2D],[5I]. Indeed, our experience in 
condensed matter and atomic physics teaches that often 
combined effects of reduced dimensionality of a quantum 
system and inter-particle interactions may effectively en- 
hances the effects of correlations and eventually may lead 



to a new, collective states of matter with properties that 
are very different from those of the underlying particles 
(one notable example is a Luttinger liquid state of inter- 
acting fermions or bosons in ID). This kind of philosophy 
is one of the driving idea of a quest for the novel corre- 
lated states of photons 22J. 
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FIG. 1: Schematic view of our proposed system to generate 
strongly-correlated states of photons: the topological pho- 
tonic insulator (bulk of the left picture) has topologically- 
nontrivial band structure (at the right picture). If the total 
Chern number of the bands below the gap is 1, the state in 
the gap is formed (red line in the right picture). This state 
corresponds to the chiral edge state of unidirectionally propa- 
gating photons (thick red line at the boundary of the photon- 
ics crystal) . We suggest to insert emitters of any internal level 
structure at the edge state (here only two-level structures are 
shown). Pumping the edge state from aside by a few-photon 
source (in-state) creates strongly-correlated photonics states 
at the outcome of the structure (out-state). 

One of these collective phenomena is a topological in- 
sulating/superconducting state. Recent progress in ex- 
perimental construction of topological insulating states 
for electrons [53] has generated a bunch of experimen- 
tal and theoretical activity in search for the novel topo- 
logical states of matter [24]. A central signature of the 
quantum state to have a topological property is an exis- 
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tence of a chiral edge state, a unidirectionally propagat- 
ing waves/particles which is insensitive to the local per- 
turbations (such as imperfections, local potentials, etc.). 
The whole existence of these states is due to non-trivial 
topological properties of the band structure in the bulk 
of a material. 

Recently a unidirectionally propagating chiral states of 
electromagnetic waves have been suggested theoretically 
in |25j and observed experimentally in [26j which gen- 
erated a subsequent large theoretical and experimental 
activities [37] , [25J , [2S] . These states are an optical ana- 
logue of the quantum Hall edge states. It is known that 
Hall effects belongs to a broader class of topological in- 
sulators exhaustively classified recently [30] . These topo- 
logical states are characterized by nontrivial Chern num- 
ber in the bulk and chiral topologically protected edge 
states. In 2D photonic topological insulators the edge 
states form a ID waveguide with unidirectional propaga- 
tion of photons, while in 3D topological photonics insula- 
tors the edge state is a 2D surface state with Dirac struc- 
ture of a spectrum [51] . Topological insulators for pho- 
tons can be constructed by engineering specific magneto- 
optical photonic crystals producing a band structure with 
nontrivial Chern numbers. 

Here we suggest to use the unidirectional edge states 
(chiral edge state in 2D topological photonic insula- 
tors) as a robust platform for controllable generation of 
strongly-correlated photonic states. To achieve this goal 
we suggest to couple an array of emitters to the edge 
state. Photons in the edge state populated by an ex- 
ternal few-photon source interact with an ensemble of 
emitters and produce an outgoing photonic states with 
robust, controllable and universal properties. An origin 
of these robustness and universalities lies in the topologi- 
cal nature of the edge state. We found in particular that 
the outgoing photonic wave function has no information 
about position of emitters while its zeroes are given by 
the universal numbers - zeroes of Laguerre polynomials. 
Even on a few-photon level the wave function has "frac- 
tionalized" form and position of one zero imply the others 
thus making the outgoing state strongly correlated. We 
also observe strong even-odd effect with respect to the 
number of emitters corresponding to switching between 
bunching and anti bunching behavior for even and odd 
number of emitters respectively. We show also that all 
these properties are robust with respect to effects coming 
from temperature, imperfections and/or disorder. These 
correlated states of photons can be used in emerging 
quantum devices. 

ID EDGE STATES WITH EMITTERS 

We consider unidirectionally propagating topologically 
protected chiral edge electromagnetic states in photonic 
crystal. We assume that an array of two- or multi-level 



systems can be inserted at that chiral edge channel. We 
assume that transition frequencies of that systems (which 
we call emitters) are close to the frequencies of the prop- 
agating chiral mode. Moreover, since the origin of the 
chiral edge comes from the Dirac-like cone structure in 
the k-space, it is reasonable to assume a linear dispersion 
for the propagating chiral model. Therefore, our model 
has a form of electromagnetic waves with a linear disper- 
sion, propagating in one-dimensional edge geometry and 
interacting with a system of two- or multi-level emitters. 
The Hamiltonian can be written therefore as 

M 1 

H = -i dxa) (x)d x a{x) + +^ A b {S§ + -) 

J b=l 

M 

- ^$> fc + a(r h ) + at(r b )S b -] (1) 

6=1 

Here, the electromagnetic field operators a){x), a(x) sat- 
isfy a usual commutation relations [a(x), ^(y)] = S(x—y) 
and have a linear dispersion spectrum, whereas the 
emitter-related variables S£ are placed at coordinates r a 
at the ID edge. Emitters are determined by the inter- 
nal level structure and by the allowed transitions between 
them. Thus, if we have a two-level system then a — x,y, z 
and = S x ± iS v ) satisfy the usual spin algebra of 
atomic operators [S z , S±] = ±S ± , [S+,S-] = 2S Z and 
we have a natural representation of the transitions in the 
two-level system with transition frequency given by Ab 
for the emitter at position r b - For the array of emitters 
we can allow the frequencies to be distributed around 
some average frequency A. On the other hand, if we 
have a three-level structure of the emitter, the definition 
of operators S a depends on the type of a scheme, which 
can be of either types, A, E, or V. Thus for the A-type 
scheme we have S+ = <?3i|3)(l| + #32|3)(2| while for the 
F-scheme S + — (?3i |3)(l|+.g2i|2)(l| and for the E-scheme 
S+ = .9 32 |3)(2| + . 9 2i|2)(l|. In all cases S_ = S\_. In this 
paper we mainly focus on the two-level case. 

The problem is further specified by the initial state. 
Here we assume a natural initial state inspired by the cur- 
rent experimental realizations of waveguides with emit- 
ters coupled to it [TS],[T7]: we assume that the edge state 
is pumped in by the external few-photon source while 
the emitters are initially in the ground state. Injected 
photons unidirectionally propagate as a wave packet at 
the edge state and interact with the ensemble of emit- 
ters. The resulting state is read-off by the detector. In 
the basis of a Fock states \in) = J dkif({ki}) TJ^ a\ |0) 
where i = 1 . . . N and /({fci}) is an envelope function 
of momenta {ki}. The evolution of this initial wave- 
function is given by the scattering matrix S(N; M; {h}), 
\out) — S(N;M;{l})\in) where the scattering matrix 
S(N; M; {I}) describes a scattering of iV-photon state 
| in) on the array of M emitters each of which may have a 
different number of levels characterized by a set of num- 
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bers {lb}, h = 2, 3, 4, . . . for b = 1 . . . M. In the case of 
two-level emitters the model ([I]) can be solved exactly by 
the Bethe Ansatz and the scattering states for any initial 
conditions can be determined exactly [33] • Basic ideas 
of this approach are summarized in the Supplement. To 
solve a more general problem with combination of emit- 
ters of different level structures we are using a combina- 
tion of Bethe ansatz [32] , [33] , [34] , [35] and a scattering 
matrix formalism [3S] (see also Supplement). A combina- 
tion of these methods provides a complimentary pieces of 
insight into how the combination of several basic ingredi- 
ents of a model (position of emitters and their level struc- 
ture, strength of light-emitter interaction, etc.) serves as 
a powerful tool for engineering correlated multiphoton 
states. 

Various results of our computations suggests the fol- 
lowing fundamental property of the model: the scattering 
matrix of the TV-photon state on the array of M atoms 
is given by the convolution of the ./V-photon scattering 
matrices on every atom, 

S tot (N; M: {l b }) = S(N; 1; k) * ... * S(N; M; l M ) (2) 

where the convolution (*) can be taken in any convenient 
representation (e.g. coordinate or momentum spaces). 
Sufficient conditions for this property to be correct are 
the following: (i) unidirectional nature of a spectrum, 
and (ii) a constant group velocity of the incoming wave 
packet. Indeed, in the case of a chiral spectrum the trans- 
fer matrix is identical with the scattering matrix provided 
that the group velocity is a constant for all particles. As 
a first straightforward consequence we conclude that the 
outgoing state of photons does not depend on the posi- 
tions of emitters. This is the origin of many universal 
properties of the outgoing photonics states. [412] 



Few-photon initial conditions 



where z = x + 1, t is time. Here 

A a - A b - in 



A„ - Aft 



(4) 



First we consider the limit Aj — > A, which can be per- 
formed analytically (see Supplement). The calculation 
yields 



Mv) = 5 (y~ z ) 



y))e" 



(5) 

- K (z-y)/2 



where L^_ 1 {x) is an associated Laguerre polynomial. 
Assuming Gaussian form of the incoming wavefunction 
0m(aj) = cr- 1 / 2 ^- 1 / 4 exp [i{A + S)x - x 2 /2a 2 } , where 
6 is the detuning w.r.t. the atomic transition fre- 
quency A, the outgoing wavefunction can be written as 

<?Wt(y) = M{y — z) + 0i,scatt(y - z). 

In Fig. pj we plot |0i jScat t| 2 for different M and de- 
tunings S. For different set of parameters the plots are 
shown in Supplement. It should be remembered that 
\(/>out(y)\ 2 = |^i, sca tt(j/)| 2 only for \y\ > a, since in the 
region \y\ < a the scattered part is masked by the non- 
scattered part. 

These results manifest clearly a strongly-correlated 
and universal character of the scattering in our system. 
The zeroes of the Laguerre polynomials determine the 
positions of minima of the outgoing wave packet. Thus 
if the position of the first peak is known the subsequent 
peaks and minima are universally determined by the ze- 
roes of Lw_ 1 (x). The result does not depend on positions 
of M emitters and thus universal. Time delay on every 
emitter formed by the emitter-photon bound state in Q 
(exponentially decaying part of the wave function) leads 
eventually to "fractionalization" of the wave packet into 
M — 1 peaks. 



We discuss scattering of a few-photon wave packets on 
arrays of M emitters. If we measure all frequencies in 
units of k (thus e.g. A p h ys — kA) and the distances in 
units of ( e.g. y v h ys — v/k), the coupling constant 
enters the problem exclusively through this renormaliza- 
tion of frequencies and distances. 



Single-photon scattering 

Using the formula derived in [34j . we investigate the 
scattering of a chiral photon on array of M atoms of tran- 
sition frequencies Ab, b = 1 . . . M. The formula for the 
out-going wavefunction produced by the S (x) -function in- 
state reads 

Mv) = S(y~z) (3) 
- K6(y<z)^2c a exp[(iA a + K/2)(y-z)]. 



Two-photon scattering 

The two-photon scattering can be shown to have three 
parts: (i) the one corresponding to no interaction with 
emitters, (ii) the one corresponding to one photon in- 
teracting with emitters while a second one goes through 
the ensemble without interaction and (iii) the irreducible 
part of scattering corresponding to two photon-bound 
state formed in the presence of emitter. This picture 
of different scattering channels will be important below 
to understand the outcome of the scattering. 

Let us look into the Gaussian incoming wavefunction 
of the form 



1 



: exp[i(A + S)(xi + x 2 )] 



x exp[-(xi + x 2 ) 2 /8fj 2 - (n - x 2 f/2a 2 ], 
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result up the first nontrivial order S 2 reads 
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FIG. 2: One photon scattering off M emitters for detun- 
ing S — (upper panel) and scattering off M = 10 emit- 
ters. Initial state has a form of Gaussian wave packet with 
dispersion a. For other values of a see the Supplement. All 
frequencies are expressed in the units of k and 5 p h ya = k5 and 

O'phys = K> 0" • 



where S is again a detuning from the atomic resonance 
frequency. We are mainly interested in the limit of wide 
pulse fi — > oo, however the relative position of the two 
photons is confined into the interval of the order of a. 
When the limit Af, — > A is taken (see Supplement) we 
observe again a polynomial structure related to Laguerre 
polynomials. Moreover, the scattering in the case of 6 = 
does not depend on the number of atoms, only of its 
parity. These limiting forms of the outgoing wavefunction 
reads 



i(d) 



for even M even, and 



1 



e 2 ^ 



— V 2-Kae 2 + 8 crfc 



2V2 



(7) 



(8) 



for odd M. Here d = yi — 1/2 is a photon's relative coor- 
dinate. One can also calculate the asymptotic expansion 
of the outgoing wavefunction for large detuning i5. The 
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FIG. 3: Two photon scattering on M-atoms. Note that the 
detuning 5 and the width of the pulse a are measured in units 
of k, i.e. Sphys = k<5 and a p h ys = k~ x o. Note, that while for 
M even the photons do not scatter at S = 0, for M odd they 
do. Note that for large detuning 8 the effect of scattering 
asymptotically vanishes. 

There are several effects present in the scattered state 
off M emitters: one is a renormalization of the incoming 
Gaussian waveform, which is likely not to be detectable, 
the other is the exponential tail, which scales as S~ 2 re- 
gardless of the number of atoms. This is another univer- 
sal feature of the scattering in our system. The exponen- 
tial tail is apparently a remnant of single-particle-only 
scattering, which can be seen either by direct compari- 
son with 1-particle scattering formulas, or by employing 
the scattering matrix in fc-space. There one can easily 
see that the irreducible two-particle part of the scattering 
scales as a higher power of 8 than the one-particle scatter- 
ing. Further, the outcome of the two-particle scattering 
is sensitive to the parity of the number of emitters: for the 
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odd number of emitters a tendency towards anti bunch- 
ing is pronounced while for the even number of emitters 
outgoing photons are bunched, see Figs The poly- 
nomial nature of the scattered wave function persists for 
the two-particle scattering as well. 

Robustness of correlated states 

Next we consider a more physical situation when 
atomic frequencies are random variables distributed ac- 
cording to the Gaussian distribution. The source of the 
randomness can be manyfold. In particular it can be 
caused by the temperature since for existing topological 
photonic crystals operating in the GHz regime the tem- 
perature effects can be important. This could also be 
imperfections or external noise coupled to emitters. We 
will model the broadening of the levels by the random 
Gaussian ensemble with a variance £, so that the distri- 
bution of levels ~ exp(- A2/2£ 2 )<5(£ a (A a - 6)) 
is taken around some mean detuning 5. Based on ex- 
act expressions we can perform the ensemble averaging 
of the observable quantities. In Figs. Q we plot typical 
results for the averaged two-photon probability density 
corresponding to correlation function. We observe ro- 
bustness of the results we described above for the big 
range of variations of parameters. We conclude that our 
results are robust with respect to level's broadening. 

DISCUSSION 

We suggested a novel approach of producing strongly- 
correlated states of photons in a highly controllable way. 
Our setup uses the edge states of photonics topological 
insulators with multi-level emitters coupled to it. We 
observed a number of universal features of scattering in 
our setup: multi-particle scattering does not depend on 
position of emitters. The outcome of the scattered wave 
packet has a polynomial structure and the minima of 
the outgoing pulse are given by zeroes of this polyno- 
mial. For a single-particle (and reducible part of multi- 
particle scattering in general) case these polynomials are 
Laguerre polynomials. In a single-particle case the scat- 
tering thus looks like if the photon would be "fraction- 
nalized" in between the minima of that polynomial. 

Physical realization. To observe quantum many-body 
effects of photons described here one needs a photonic 
topological insulators introduced and studied in [35] , [27J . 
The role of emitters can be played by the quantum dots or 
superconducting qubits, like e.g. in [TO] and [J5]. Also, 
available single-photon emitters made of quantum dots 
coupled to photonic crystal waveguides makes our system 
realizable at the today's level of development [TT] . [15] . [5] . 
Typical frequency ranges of existing photonic topologi- 
cal insulators are GHz with matches the photonics single 
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FIG. 4: Averaged two-photon scattering off M emitters. The 
transition frequencies A a are independent random variables 
distributed from the Gaussian distribution of variance E. The 
detuning <5 is the deviation of the incoming wave frequency 
from the sample mean frequency. The solid lines correspond 
to the ensemble average of the density (second-order correla- 
tion function). The filled regions denote the median (darker) 
and mean (lighter) absolute deviation. The detuning 5, the 
width of the pulse a and the variance E are measured in units 
of k, i.e.Sphys = k5 and a p hys — k~ 1 <j. 

emitters frequencies and couplings. Temperature effects 
relevant for this frequency range show up in the inhomo- 
geneous broadening of levels. We take this into account 
by assuming Gaussian fluctuations of detunings for ev- 
ery emitter. Photon's initial state injection can be made 
using the same single-photon emitters. 

Extension for a general level's scheme. Here we fo- 
cused on a two-level structures, however using the con- 
volution property[2]one can immediately apply several re- 
cent results on scattering off three- or four-level schemes 
[201 |21j . which contain a number of interesting experi- 
mentally observable properties [37J including singlc- 
cmittcr-EIT, transistor-like behavior, implementation of 
quantum gates (e.g. CNOT [38]), etc. By combining 
emitters of different level structure this flexibility makes 
it possible to build an optical schemes with desired cor- 
related properties of outgoing photonic wave function. 
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For explicit evaluation of outgoing state one should use a 
scattering formalism [36] (see also Supplement) together 
with the convolution property ph. 

Extension for different initial conditions. Here we ex- 
plicitly evaluated the properties of the scattered light for 
photonic initial conditions, when emitters are not excited 
initially. Other initial conditions can be treated similarly. 
In particular one can think of emitters initially excited 
from aside. In this case the evolution can be studied us- 
ing similar method while the outgoing wave function will 
depend on the emitter's positions. Also coherent state 
initial condition can be analyzed using our experience 
with one- and two-photon scattering as well as a knowl- 
edge of scattering on a single emitter. If we focus on the 
case of a single- mode fco coherent state \a)k a than it be- 
comes clear that a multi particle irreducible parts make a 
minor contribution to the outcome \out)k in momentum 
state kg. One can draw this conclusion by expanding 
the coherent state into a Fock basis and using the ob- 
servation that multi-particle irreducible part has smaller 
available phase space than the reducible parts. Observed 
universalities of a single-particle scattering will be there- 
fore present in the coherent state initial state. More de- 
tailed arguments are given in the Supplement. 

Effects of dissipation of photons into outside world is 
equivalent to opening of additional scattering channels in 
our approach. Although detailed studies of decoherence 
of topologically-protected states are not available to our 
knowledge, based on experimental observation of long- 
lived edge states , [37] one can generically believe that 
topologically protected states will be robust with respect 
to dissipation. 

Currently available topological photonic crystals op- 
erate at GHz frequency range. This is related to diffi- 
culties of creating large magneto-optical response (which 
is needed for nontrivial topology of the band structure) 
with existing materials. However some recent develop- 
ments [39 in material science may lead to significant 
magneto-optical response in the optical domain which 
then may be used to engineer topological photonic crys- 
tals in the optical range of frequencies. 

Finally we note that in 2D interacting photons can lead 
to anyonic statistics and fractional Hall states [40]. Re- 
cent suggestion of realizing Dirac cone structure [3T] may 
pave the way to this exotic physics of interacting 2D pho- 
tonics fluids. Moreover recent classification of interact- 
ing bosonic topological insulators [U] suggest a possible 
route to engineer even more exotic states of light. 

We expect that the novel correlated states of photons 
or plasmons will find a number of applications, some of 
them we can only envision now: quantum information 
technology, optomechanics and precision measurements. 

Acknowledgement. We would like to thank V. Yudson 
for useful discussions. M.R. and V.G. are supported by 
MaNEP and Swiss NSF. M. P. is supported by DFG-FG 
723. V.G. thanks KITP for hospitality. 



[1] K. Vahala, (ed.) "Optical Microcavities" , (World Scien- 
tific, Singapore, 2004). 

[2] M. Fleischhauer, A. Imamoglu, and J. P. Marangos, 
"Electromagnetically Induced Transparency: Optics in 
Coherent Media", Rev. Mod. Phys. 77, 633 (2005). 

[3] K. Hammerer, A.S. Sorensen, E.S. Polzik, Rev. Mod. 
Phys. 82, 1041 (2010) 

[4] I. Carusotto, C. Ciuti, |arXiv: 1205. 6500| (2012). 

[5] G. Wrigge, et al., Nature Phys. 4, 60 (2008). 

[6] J. Hwang et al., Nature 460, 76 (2009). 

[7] A. A. Abdumalikov, et al. Phys. Rev. Lett. 104, 193601 

(2010) . 

[8] M. D. Eisaman, J. Fan, A. Migdall, and S. V. Polyakov, 

Rev. Sci. Instrum. 82, 071101 (2011). 
[9] T. Volz, et al. Nature Photonics 6, 605 (2012). 

[10] C. Lang, et al. Phys. Rev. Lett. 106, 243601 (2011). 

[11] P. Yao, V. S. C. Manga Rao and S. Hughes, Laser Pho- 
tonics Rev. 4, 499 (2010). 

[12] K. P. Nayak, K. Hakuta, New J. Phys. 10, 053003 (2008); 
K. P. Nayak, Fam Le Kien, M. Morinaga, and K. Hakuta, 
Phys. Rev. A 79, 021801(R) (2009); G. Brambilla, J. 
Opt. 12 043001 (2010); Fam Le Kien, S. Dutta Gupta, 
K. P. Nayak, and K. Hakuta, Phys. Rev. A 72, 063815 
(2005); G.Sague, E. Vetsch, W. Alt, D. Meschede, and 
A. Rauschenbeutel, Phys. Rev. Lett. 99, 163602 (2007). 
Fam Le Kien, V. I. Balykin, and K. Hakuta, Phys. Rev. 
A 70, 063403 (2004); V. I. Balykin, K. Hakuta, Fam 
Le Kien, J. Q. Liang, and M. Morinaga, Phys. Rev. A 
70, 01 1401 (R) (2004). K. P. Nayak, P. N. Melentiev, 
M. Morinaga, Fam Le Kien, V. I. Balykin, K. Hakuta, 
arxiv:quant-ph/0610136 (2006); I. Schuster, et al. Nature 
Physics 4, 382-385 (2008); Dayan, B. et al. Science 319, 
1062-1065 (2008); Fink, J. M. et al. Nature 454, 315- 
318 (2008). Hoi, I.-C. et al. Phys. Rev. Lett. 107, 073601 

(2011) . 

[13] A. Wallraff, et al. Nature 431, 162 (2004). 

[14] T. Aoki, et al. Nature 443, 671 (2006). 

[15] O. Astafiev, et al. Science 327, 840 (2010). 

[16] J. Claudon, et al. Nature Photonics 4, 174 (2010); 

[17] A. Reinhard, et al. Nature Photonics 6, 93 (2012). 

[18] I-C. Hoi, Phys. Rev. L ett. 107, 073601 (2011); I.-C. Hoi, 
et al. |arXiv: 1210.4303 (2012). 

[19] H. F. Hofmann, K. Kojima, S. Takeuchi, K. Sasaki, Phys. 
Rev. A 68, 013803 (2003); J. T. Shen, S. Fan, Opt. Lett. 
30, 2001 (2005); J.-T. Shen and S. Fan, Phys. Rev. Lett. 
95, 213001 (2005); Phys. Rev. A 76, 062709 (2007); D. E. 
Chang, A. S. Sorensen, E. A. Demler, M. D. Lukin, Na- 
ture Physics 3, 807 (2007); J-Q. Liao, C. K. Law, Phys. 
Rev. A 82, 053836 (2010); T. Shi, S. Fan, and C. P. Sun 
Phys. Rev. A 84, 063803 (2011); 

[20] T. S. Tsoi and C. K. Law, Phys. Rev. A 80, 033823 
(2009); D.Witthaut and A. S. Sorensen, New J. Phys. 12, 
043052 (2010). H. Zheng, D. J. Gauthier, H. U. Baranger, 
Phys. Rev. A 85, 043832 (2012). 

[21] D. Roy, Phys. Rev. Lett. 106, 053601 (2011); D. Roy, 
Phys. Rev. B 81,155117 (2010). 

[22] M. J. Hartmann, F. G. S. L. Brandao, M. B. Plenio, Na- 
ture Physics 2, 849 (2006); A. D. Greentree, C. Tahan, 
J. H. Cole, L. C. L. Hollenberg, Nature Physics, 2, 
856 (2006); D. E. Chang et al.Nature Physics 4, 884 
(2008); I. Carusotto, et al. Phys. Rev. Lett. 103, 033601 



7 



(2009); A. L.C. Hayward, A. M. Martin, A. D. Greentree, 
Phys. Rev. Lett. 108 , 223602 (2012); C.-E. Bardyn, A. 
Imamoglu, arXiv: 1204. 1238 (2012). 



[23] M. Konig, et al. Science 318, 766 (2007). 

[24] M. Konig, et al. J. Phys. Soc. Japan 77, 031007 

(2008) ; X.-L. Qi, S.-C. Zhang, P hysics Today, Jan 2010, 
p. 33; M. Z. Hasan, C. L. Kane, |arXiv: 1002.3895| J. E. 
Moore, Nature 4 64, 194 (2010); X.-L. Qi, S.-C. Zhang, 
larXiv:10 08.2026, 

[25] F. D. M. Haldane, S. Raghu, Phys. Rev. Lett. 100, 
013904 (2008); S. Raghu, F. D. M. Haldane, Phys. Rev. 
A 78, 033834 (2008). 

[26] Z. Wang, Y. Chong, J. D. Joannopoulos, and M. Soljacic, 
Phys. Rev. Lett. 100, 013905 (2008); Nature 461, 772 

(2009) . 

[27] S. R. Zandbergen, M. J. A. de Dood, Phys. Rev. Lett. 
104, 043903 (2010); C. He, et al. J. Appl. Phys. 107, 
123117 (2010); J. -X. Fu, R. Liu, Z. Y. Li, Appl. Phys. 
Lett. 97, 041112 (2010); Y. Poo, et al. Phys. Rev. Lett. 
106, 093903 (2011); J.-X. Fu, et al. Phys. Rev. Lett. 98, 
211104 (2011). 

[28] Z. Yu, G. Veronis, Z. Wang, and S. Fan, Phys. Rev. Lett. 

100, 023902 (2008); K. Fang, Z. Yu and S.Fan, Phys. 

Rev. B 84, 075477 (2011). 
[29] T. Ochiai, Journal of Physics: Cond. Matt., 22, 225502 

(2010) . 

[30] A. P. Schnyder, S. Ryu, A. Furusaki, A. W. W. Ludwig, 
Phys. Rev. B 78, 195125 (2008); AIP Conf. Proc. 1134, 
10 (2009); New J. Phys. 12, 065010 (2010). 

[31] J. Bravo-Abad, J. D. Joannopoulos, M. Soljacic, 
|arXiv:1204.0500| (2012); L. Lu, L. Fu, J. D. Joannopou- 



los, M. Soljacic, [arXiv : 1207.0478 (2012). 
[32] V. I. Rupasov and V. I. Yudson, Sov. Phys. JETP 59, 
478 (1984). 

[33] V. I. Yudson, Sov. Phys. JETP 61, 1043 (1985). 
[34] V. I. Yudson, Phys. Lett A 129, 17 (1988). 
[35] V. I. Yudson, P. Reineker, Phys. Rev. A 78, 052713 
(2008). 

[36] M. Pletyukhov and V. Gritsev, New. J. Phys. 14, 095028 
(2012). 

[37] L. Slodicka, G. Hetet, S. Gerber, M. Hennrich, R. Blatt, 

Phys. Rev. Lett. 105, 153604 (2010); Miicke, M. et al. 

Nature 465, 755-758 (2010). T. Kampschulte, et al. Phys. 

Rev. Lett. 105, 153603 (2010). 
[38] H. Zheng, D J. Gauthier, H. U. Baranger, 

|arXiv:1211.l7Ti"1 (2012). 
[39] V. I. Belotelov, et al. Nature Nanotechnology 6, 370 

(2011). 

[40] R.Y. Chiao, T.H. Hansson, J.M.Leinaas, S. Viefers Phys. 

Rev. A 69, 063816 (2004). 
[41] X. Chen, Z.-C. Gu, Z.-X. Liu, X.-G. Wen, Science 338, 

1604 (2012). 

[42] Note however, that one has to distinguish two different 
physical situations. If the typical wavelength of the in- 
coming state is larger than the average distance between 
some group of m < M atoms, a Dicke-like collective m- 
emiter is formed. In this case for the two-level systems 
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Supplementary Information 



EXACT SCATTERING STATES 



Yudson's solution 



A correct basis of the scattering states for the model introduced in Eq. ([I]) is constructed using polaritonic operators 

M „ t 
b=l b 

where (6=1... M) are coordinates of emitters and numbers A — > Xj (j — 1, . . . , N) characterize the excitations in 
a system of photons coupled to M emitters. A computation of the time evolution of the initial states can be performed 
using the method of Bethe ansatz developed for this and similar problems in |33j . Using this resummation method, 
the evolution of the initial A^-photon state of the form |^q) = Tlj=i a ( x j)\fy i s obtained as an exact solution of the 
evolution problem. In the asymptotic limit of large time only photonics components of polaritons ([9| remain and 
therefore the outgoing state has the form 

N 

\out)= I d N y^ N (y)l[a^y j )\0) (10) 



where 



d> UA I dNX TT A '~ A J+^[ si S n (^~^)] (-.is 




Here = scj + i and the integration contour T is determined as follows: in the complex plane of parameters A the 
contour is a composition r = 71 © 72 ... © "/n such that: (i) hnrfj+i-bar/j > k for j = 1, . . . , N and (ii) at Re7j = A a 
one has Inryj > —k/2 for Xj < r a . Here we assume that initially photons are injected in the most left region of space 
with respect to all emitters and that all emitters are initially in their ground states. One can also consider more general 
initial conditions, when some (or all) emitters are initially in the excited states. In this case the condition (ii) for the 
choice of the contours jj should be modified by adding two more: Imr/j < k/2 for Xj > r a and — k/2 < Imr/j < k/2 
for Xj = r a . Note that here we corrected the definition of contour T with respect to Ref. [53] . 
Let us consider several important examples. 



Single-photon scattering 

First we consider a scattering of a single chiral photon state off an array of two-level systems. The evolution of the 
plane wave state \k) = J dxe lkx E + (x)|0) = E + (k)\0) can be computed using the general scheme of the evolution of 
the basis state. In this case the amplitude 4>\{y) can be computed 34J as (this result is also quoted in the main text) 

My) = S$ = S(y - z) - K6{y < z) (13) 



^C a exp[(iA Q + K /2)(y-z)] (14) 



where 



Ca= n A Q - Ah • (i5) 

The out-state of the wave is therefore given by 

\k)out =t(k)E+{k)\0) (16) 
where (repeating the steps of section III. A of [35]) 

Tjk-Aa-in/2 

f w = n Jb _ Aa+<lc/2 ' 

a ' 

which reproduces known results for the chiral part of scattering state [19]. It is a convolution (product) of scattering 
phases on individual emitters |36j . 

In the special case of equal detunings A a = A the scattering phase (17) turns into (fc — A — ik/2) m / (k — A + ik/2) m , 
and therefore 

My) - i I ( 1 - ^h^Y - S(y -.) + £ c;:t^l f ^ 



2tt J V k-A + in/2 vtf ' ^ m 2n J (fc — A + in/2) 

m—1 ' ' ' 



5(y -z)- ,9(z - y) f) (% ti&^^ e «* + */m-.) 



m—1 

-(i) 



= S(y- z)- KQ(z-y)L^ ) _ 1 (K(z~y))e {lA+K/2 '> {y ~ z \ (18) 
where L^_ 1 (x) are associated Laguerre polynomials. 
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FIG. 5: One photon scattering on M-atoms for detuning S — and the initial Gaussian wave packet (see the main thext) 
characterized by a. 



Two-photon scattering 

(2) 

By definition, the two-particle scattering matrix S M scatters the in-state 

| in) = / dz l dz 2 tp(z l ,z 2 )a 1i (z l )a jl (z 2 )\0), 



normalized as 2 J dz\dz 2 \il)[z\, Z2)\ 2 = 1, into the out-state 

I out) = J dy 1 dy 2 dz 1 dz 2 ip(z 1 ,z 2 )^S ( ^ (y 1} y 2 ; z 1} z 2 )a? (y^a^ (y 2 )\0) . 
For the case of scattering off an array of M emitters the scattering matrix reads 



(19) 



(20) 



f 2) 



d\i f d\ 2 



71 2n 7 72 2n 



1 - 



M 

n 



Ai — A a — in/2 



M 

n 



\ AJ : Ai - A a +w/2. , 

\a=l ' / \6=1 

+ (2/1 <-> Vi) ■ (zi <-» z 2 ) 



2inQ(y 2 - t/i) 
Ai - A 2 + ik 

X 2 — — in/2 
A 2 - A b + in/2 



^l{vi-zi)+i\ 2 (y2-z 2 ) 



SM(yi,zi)S$(y2,z 2 ) + S ( g(y 1 ,z 2 )S ( g(y2,z 1 )+iT^\yi,y2,zi,Z2), 



(21) 



where Sj9 (x, y) are the single-particle scattering matrices. The first two terms in the last line represent the reducible 

(2\ 

part of the scattering (independent single-particle scattering) while TJ^ (2/1,2/2, -Zi, -^-matrix represent correlated 
irreducible part. The reducible part of the two-photon scattering off M emitters can be written as 



S^U^Ay^z) 



1 

2^ 



dpdke ipAy - lkAz [S(p -k) + 5(p + k)] 



5 + k-in/2\ M (5-k-in/2 K Ai 



5 + k + in/2J \S -k + in/2 



In case of resonance S = we obtain 



SM;red(°> A 2/> Az) = ^ / dpdke lpAy ~ lkAz [S(p - k) + S(p + k)] = S(Ay - Az) + S(Ay + Az). 



(22) 



(23) 



It is also convenient to define 



$M £d(0; Ay, Az) = Q(Ay)Q(Az)S$M-, Ay, Az) = 6(Ay - Az) 



(24) 
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The irreducible two-photon part of the scattering is given by 



(2) 

Tm {Vi, V2,zi, Z2) = -2kQ(zi - Z2)Q {V2 - yx) 



M 

n 



71 27T J l2 27T 

M 



Ai — A2 + in 



Ai — A a — in/2 ^ j j-j A2 — Ah — in/2 



iAi(yi -21 )+iA 2 (2/2-22) 



\ 1 \ Ai - A a + in/2 / \i A2 — A;, + in/2 

\a=l ' / \b=l / 

+ (yi 2/2) ■ (zi 22)- 



(25) 



Note that for the calculation of this quantity the mutual arrangement of the contours 71 and 72 does matter. Following 
the above prescription for a choice of integration contours one can show that (cf. Eq. (21) from Ref. |34| ) 



T^\yi,y2,z 1 ,z 2 ) = -2k 3 6(zi > z 2 > y 2 > J/i)^ 



C„C\ 



a,b 



A n . - Ah + in 



e i(A a -iK/2)(y 1 -z 1 )+i(A b -iK/2)(y 2 -z 2 ) 



+{yi <"> 2/2) ' (Zl 4-> Z 2 ). 



(26) 



A special case of M equal detunings is interesting. Let us consider the limit A a ^ — > A of (26) for Z\ > z<x > 2/2 > y\- 
At first, we introduce an integral representation 



T$ = 2ik 3 dre~ KT ^C„e« A »-<" ■W(r+ux-i) S j (^C b e 



i(A b -iK,/2)(-T+y 2 -z 2 ) 



which helps to decouple the sums over a and b. At second, we employ the identity [cf. Eq. (j 1 8|) ] 

M 

I dse.- tsx 

2tt 



i(A a -iK/2)x 



M 

— / dse~ lsx V C M - 
2ir J ^ m ( 



(-«)' 



With its help Eq. ( 26 ) is transformed in this case into 
- (2 ) ~. f ds f dr f Zl - yi 



2tt J 2tt J 

M M 

X C mi C : 
mi — 1 m 2 — 1 

ds f dr 
2tt J 2tt 

\ i M 

X 

m\ — l m 2 =l 



^ (s - A + ^/^/2) r 



^ Te -i(r-s-iK)r e -is(2i-t/i) e -ir(2 2 -j/ 2 ) 



(27) 



(28) 



|'»2 



2lK 



A/ H*)" 1 

m2 (s- A + m/2)™i (r - A + in/2) m 2 

^ Te -i(s+r-iri)T e -i(s+iK)(yi-zi) e -ir{z 2 -y 2 ) 



EST r M c M 
/ ; "I "12 



(«) r 



{-in)' 



(s + A + m/2) m ^ (r - A + inj2) m ^ 



Integrating over r we obtain 



r 1 



■(2) 



A/ 



2 k 



2 k 



M M 



r — s — IK 



2tt 7 2tt 

(is /" <^ r e -K(^l-I/l) e -*I , (zi-l/l+22-!/2) 

s + r — ir/ 



^2 X! ^mi^i 



(-<«) 



m 2 



mi = l 7712 = 1 

M M 



AT (~^) T 

™ 2 (s - A + in/2) m i (r - A + «k/2)" 1 2 



(29) 



(30) 



v v r M r M (^) mi H«) m3 

2^ mi ™ 2 (s + A + iK/2) m i (r — A + i/s/2)"*" 

mi — 1 mo — 1 



2?r / 2?r 



This expression is remarkable, since it resembles Eq. (26). However, instead of the summation over distinct discrete 
detunings A a and Ah, we have now the double integral over continuous variables s and r. The second term exactly 
compensates the contribution coming from the pole s — r — in in the first term, so one can rewrite (30) as 

(-m) mi 



r (2) 

1 M 



-2k 



ds f dr e-^i-^e-^-^) 



M AI 



71 2lT K 27r 



s — r + ik 



EST c M c M 



(-«)' 



m-i —1 mo — 1 



(s - A + iK/2) m i (r - A + i«/2) m 2 



. (31) 



where the contours 71 and 72 are small circles clockwise embracing the points s = A — ifc/2 and r = A — iz-c/2, 
respectively (cf. Eq. (44) and (45) of Ref. [33]). 
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Restoring a complete form of Eq. (31) by adding the terms (y% ■<-» j/2) ■ (^1 ^ ^2), we transform it to a mixed 
representation of the full momentum pi +P2 — k\ + = E = 2(A + 8) and relative spatial coordinates Ay = yx — Vi, 

Az = Z\ — Z2 



T£\S;Ay,Az)= £ [8; Ay, Az), 

a,/3=± 



(32) 



where the component TJ^ ++ has Ay > and Az > 0, and the other ones are obtained from it by flipping the signs 
of Ay and Az. Quantatively T^ ++ (8; Ay, Az) is 



T£ )++ (5;Ay,Az) = -2iK 



ds f dr 



8-iK,/2)(Ay+Az) 



M M 



E E c&c 



v . 27T /„/ 27T (s — r + in) (28 — s — r + in) ^— ' 

'1 '2 v yv ' m 1 = lm 2 = l 



where both contours "/[ and 7 2 embrace the origin clockwise. Noting that 



(33) 



M M 



E E 

mi — 1 7712 — 1 



M r<M ( * K ) mi ( i K Y 



(s — in) 

„M 



M 



- 1 



(r - in) M 



(34) 



we cast (33 1 to 



T^ )++ (S;Ay,Az) 



2i Ke (i5-K/2)(t\y+Az) Q2M-2 

[(M- 1)!] 2 a s M-i 9r M-i 



From the identity 
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M-l 



Q r M—l 



(r — ik) m 



s — r + in 



(M-l)! 



e - ls (Ay+Az)( s _ ^)M (r _ - k) M 

(s — r + in) (28 — s — r + in) 



,M 



s—r—0 



it follows that 

(r-i K ) M 
s — r + in 



(s + in) M 



M-l 



-(r — in) 



M-l 



(35) 



(36) 



■ Y l s\r-iK) M - 1 - l = - E S i (r- i «) M - 1 - i - S M ^ S i (r-z«)- 1 ^ 



„M 



Z=0 



/=0 



= -r™^ + P M _ 2 (r) + — , 

s — r + in 

where Pm-2{t) is a polynomial of the degree M — 2. One can then show that 



d 



M-i r 



Q r M-l 



(r — «k) 



M 



(s — r + in)(2S — s — r + i/t) 



_ (M-l)! 
r=0 _ 2(,s - 8) 



(s - 28) 



M 



■;M 



(s-28-in) M (s + in) 



M 



(37) 



(38) 



Inserting it into (35), we obtain 

2iKe (i6- K /2)X QM-1 



nr^(^x) 

where X = Ay + Az. 



(M-l)! ds M ~ 



-isX 



(s - m) M 
2(s - 8) 



(s - 28) 



M 



(s - 28 - in) M (s + in) 



\M 



s=0 



The formula (39) is capable to account for the parity effect in the case 8 = 0. We have 



T M 2)++ (8 = 0;X) 



(M-l)! ds™- 1 



-isX M-l! 

3 2 



M 



1 - 



(s — in) 
(s + in) M 



= 2ine 



.•„.-kX/ 2 1-(-1) M 



(39) 



(40) 



Combined with the result for the irreducible part this leads to explicit formulas for the parity effect described in the 
main text. 



At finite 8 we represent (39 1 as 



T^ )++ (8;X) = 2iKe {l5 - K/2)x F{S/K,KX). 



(41) 
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The function F(6/k, kX) is defined by 



5M-1 



F (d/K,KX) = j^l[e- isKX m) l 



where the function 



/(*) = 



1 



(a-i) 



M 



(M — 1)! 2(s — 5/k) 



(s-26/k) m ____ 
(s -25/k- i) M ~ (s + i) M 



(42) 



(43) 



can be replaced by the polynomial because of condition s = after differentiation in (42). One then finds that 
F(S/k, kX) is a polynomial of degree Al — 1 in kX 



M—l 



M-l 



F(6/ K , K X) = £ /( ", (0) ^ [e-^U = £ ^ M - X / W (0)(-^) 



M-2-1 



(=0 



(44) 



where C*^" 1 are binomial coefficients. We therefore conclude that the outgoing wave function for the two-photon 
case develops a polynomial form as well. 

Plots for the probability density of the two-particle scattering are shown in the main text. Here we supplement 
them for different set of parameters. 




2 3 

«(yi ) 



1.0 1.5 2.0 



3.0 



FIG. 6: Two photon scattering on M-atoms. Note that the detuning S and the width of the pulse a axe measured in units 
of k. i.e. Sphys = k8 and a p h V s = k~ 1 5. For zero detuning 5 — the contributions from the one-particle and the two-particle 
scattering cancel each other and as a result the wave goes through unscattered. Note that for large detuning 8 the effect of 
scattering asymptotically vanishes. 



SCATTERING FORMALISM 



Here we describe an application of the scattering formalism developed in [36] for computing the nontrivial elements 
of the scattering S and T matrices for iV-photonic states on arrays of two- and tree-level emitters in various positions. 
The advantage of this approach is that it has no limitations of the Bethe ansatz. In particular, one can consider an 
array of emitters of different level structure (two-level, three- and four-levels) for different emitters and for arbitrary 
emitter-dependent coupling strengths n a . The approach agrees with the Bethe ansatz in the limit n a — > k for the 
two-level scheme. 



Single-photon scattering. A generalization of (17) for the scattering of a single photon off M emitters for different 



coupling constant is obtained by multiplying corresponding phase shifts on individual emitter, 

J - J - k - A a + iK a /2 



An agreement in the limit n a — ► k is obvious. Single photon scattering matrices off 3-level emitters are listed in [3"o] 
and agree with known results |20j . 
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FIG. 7: Two photon scattering on M-atoms. Note that the detuning 8 and the width of the pulse a are measured in units 
of k, i.e. Sphys = k8 and a p hys = « _1 (7. For zero detuning 8 = the contributions from the one-particle and the two-particle 
scattering cancel each other and as a result the wave goes through unscattered. Note that for large detuning 8 the effect of 
scattering asymptotically vanishes. 



Two-photon scattering. The scattering matrix (21 1 off AI emitters is a convolution of M scattering matrices off 
individual emitters, as suggested in [3B]. Let us illustrate this on example of the two emitters M = 2. 

For different k\ and k 2 the solution is presented in [36] . Here we focus on k\ = k 2 = k. First we take the Fourier 
transform of the Bethe ansatz solution ( 26 ) 

1 



IK 

n 



dy 1 dy 2 dz 1 dz 2 e- i ^- i ^+ ik ^+ ik ^T^ ) (yuy 2 ,z 1 ,Z2) 
C a Cb 1 



(27T)2 
3 

s(pi +p2-ki- fc 2 ) 



A a - A h + in pi + p-2 - A a - A b + in 



1 



1 



1 



1 



x pi — A a + in/2 p2 — A a + in/2 J \ki — A a + in/2 k 2 — A a + in/2 
and evaluate it for the case of M — 2. We have 



C x = 



Ax 



Ar 



G 2 = 



A2 — Ai — in 
A 2 - Ai 



Omitting the factor — 6(pi + p 2 — fci — k 2 ) we evaluate the sum in (46) 



k 2 = - A Al A A , 2 . {Ct{E - 2Aj + i«)M 1 (p 1 )M 1 ( P2 )M 1 (A :i )Af 1 (fc 2 ) 
h — Ai — A 2 + ik 

+C 2 (E - 2A 2 + iK)M 2 (p 1 )M 2 (p 2 )M 2 (fei)M 2 (fc 2 )} 

i/i 

_£/ — Qfl — cv 2 / 

x [(E - 2a 1 )M 1 (pi)M 1 (p 2 )M 1 (A!i)M 1 (* 2 ) + (E - 2a 2 )M 2 {p 1 )M 2 {p 2 )M 2 {k 1 )M 2 {k 2 )] 



= K Z \l- 



IK 

Ai- A 2 



1 - 



— oi\ — a 2 



x [(E - 2a 1 )M 1 ( Pl )M 1 (p 2 )M 1 (k 1 )M 1 (k 2 ) - (E - 2a 2 )M 2 {p 1 )M 2 {p 2 )M 2 {k 1 )M 2 {k 2 )] 



(46) 



(47) 



(48) 



where Mi/ 2 (p) = l/(p— cvi/ 2 ), evi/ 2 = Ai/ 2 — ik/2. Note that the expression (48) is invariant under the exchange of 
Ai and A 2 (which also implies a\ «-> a 2 and Mi <H> M 2 ). 
Using the identities 



Mxipx) = [l + Mi(p x )(Ai- A a )]M 2 (pi), 

Mi(p2) = [l + Mi( P2 )(Ai- A 2 )]M 2 (p 2 ), 

M 2 (fct) = [l-M 2 (fe 1 )(A 1 -A 2 )]Mi(A;i), 

M 2 (fc 2 ) = [l-M 2 (fc 2 )(Ai-A 2 )]Afi(A :2 ), 



(49) 
(50) 
(51) 
(52) 
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and 



M 1 [p 1 )M 1 {p 2 ) = [1 + (Ai - A 2 )(E - Ai - A 2 + iK)M 1 {p 1 )M 1 (p 2 )]M 2 (p 1 )M 2 (j) 2 ), 
M 2 {k x )M 2 {k 2 ) = [1 - (Ai - A 2 )(E - A x - A 2 + iK)M 2 {k 1 )M 2 (k 2 )]M 1 (k 1 )M 1 (k 2 ), 



we cast (48 1 to 



k 2 (E - 2a 1 )M 1 { Pl )M 1 (p 2 )M 1 (k 1 )M 1 (k 2 ) + k 2 (E - 2a 2 )M 2 ( Pl )M 2 (p 2 )M 2 (k 1 )M 2 (k 2 ) 



2iK 3 
E — ol\ — a 2 

E- Ax- A 2 

Ai-A 2 
E- Ai - A 2 



ff-Ax^-A, 
£:-«!- a 2 (Ax - A 2 ) 

Ai — A 2 



-in" 1 + 



M 2 (p 1 )M 2 (p 2 )M 1 (k 1 )M 1 (k 2 ) 
(Ax - A 2 )(£ - 2a 1 )M 1 ( Pl )M 1 (p 2 )M 1 {k 1 )M 1 {k 2 )M 2 (p 1 )M 2 (p 2 ) 



1 
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2k 4 



M 2 (p 1 )M 2 (p 2 )M 1 (A;i)M 1 (fc 2 ) 



E — ot\ — a 2 
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Furthermore, we notice that 



Sl(*i)£i(*2) = l-in(E-2A 1 )M 1 (k 1 )M 1 {k 2 ), 
S 2 ( P i)S 2 (p 2 ) = l-iK{E-2A 2 )M 2 ( Pl )M 2 (p 2 ), 



and conclude that 



T^L 2 {Pi,P2,k 1 ,k 2 ) = S pl+P2M+k2 - {k 2 {E - 2a 1 )M 1 ( Pl )M 1 (p 2 )M 1 (k 1 )M 1 (k 2 )S 2 (p 1 )S 2 (p 2 ) 
+k 2 {E - 2a 2 )M 2 {p 1 )M 2 (p 2 )M 2 {k 1 )M 2 {k 2 )S l {k 1 )S 1 {k 2 ) 



2n A 
E — ot\ — a 2 



M 2 (jp 1 )M 2 (p 2 )M 1 (ki)M 1 {k 2 ) \ , 



(53) 
(54) 



(55) 



(56) 
(57) 



(58) 



which agrees with the result of convolution [3S] evaluated at K\ = K 2 . This illustrate the convolution property ^ 
from the main text. 



SCATTERING OF COHERENT LIGHT 



Recently we studied the scattering of the initial coherent state off a single emitter [36]. Our results suggest the 
following picture of many-body effects in the coherent state. The probability of finding n& particles in the mode k 
before scattering is given by the coherent state Poisson distribution. 

P(n k ) = ^.e"", (59) 
n k \ 

where n = \a k \ 2 . The initial state reads 

|m) = ^2f3(m)\m € k) (60) 

m 

and |/?(?ti)| 2 = P(m). After scattering the state reads 



\out) = S\in) = f3(m)S\m € k) = (3(m) s(to; l)\m; m — I G k), 

m m>0 m>l>0 



(61) 
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where coefficients s are hereby defined in such a way that the state |m; I G k) with Z-particles in mode k are normalized 
to 1. 

In fact, the particles leave the mode k only if they scatter in an i-particle irreducible way with i > 2, since only then 
the photons can redistribute the total momentum between themselves and the chance to stay in k is measure zero. 
A sole 1-particle scattering conserves momentum of a photon and hence the photon cannot leave the mode k. The 
coefficients s(m; I) therefore measure amplitude of relative weight of 1-particle scattering to more particle irreducible 
scattering. E.g. s(m; 0) is generated by purely 1-p scattering (1-PIS), s(m; 1) = because of conservation of momenta, 
s(m; 2) is contributed only by 2-PIS, while s(m; 4) is contributed by both 4-PIS and combination of two 2-PIS. 

Hence we can qualitatively understand the behavior of the probability distribution P(n k ) after scattering. Indeed 



For large 7 the more-PIS dominate over 1-PIS and s(m; I) is small for small I. Morover, the states with low n k 
acquire contribution from all higher lying states. Therefore, as a result, the weight shifts to low n k which implies that 
few-particle effects are mostly important. 




(62) 



l>0 



